Optical conductivity in the normal state fullerene superconductors 
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We calculate the optical conductivity, o~{uj), in the normal state fullerene superconductors by self- 
consistently including the impurity scatterings, the electron-phonon and electron-electron Coulomb 
interactions. The finite bandwidth of the fullerenes is explicitely considered, and the vertex corection 
is included a la Nambu in calculating the renormalized Green's function. is obtained by 

calculating the current-current correlation function with the renormalized Green's function in the 
Matsubara frequency and then performing analytic continuation to the real frequency at finite 
temperature. The Drude weight in a(u>) is strongly suppressed due to the interactions and transfered 
to the mid-infrared region around and above 0.06 eV which is somewhat less pronounced and much 
broader compared with the expermental observation by DeGiorgi et al. 
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I. INTRODUCTION 

The optical spectra of fullerene superconductors in the 
normal state were found to exhibit some unusual features 
The optical conductivity, cr(u>), deviates consider- 
ably from the simple Drude behavior expected for con- 
ventional metals: the spectral weight of the Drude peak 
is reduced by about an order of magnitude and trans- 
fered to a mid-infrared (MIR) region around 0.06 eV. 
This suggests that the strong interaction effects due to 
the Coulomb and electron-phonon interactions should be 
important in the optical spectra of the fullerenes. Under- 
standing this unusual behavior in the optical conductivty, 
therefore, could reveal important information about the 
fullerenes and contribute to understanding other physical 
properties of the material. 

The optical conductivity <j(oj) represents the rate at 
which electrons absorb the incident photons at energy 
w, and is a useful probe in determining electronic char- 
acteristics of the material under study. For an ideal 
free electron gas, where the interactions between the 
electrons, and between the electrons and phonons are 
neglected, and the impurity scattering rate 1/r — » 0, 
a(ui) collapses to a delta-function, a(ui) = D tot 5(ui), 
where the coefficient D to t represents the total spectral 
weight. In this case, the optical conductivity sum rule, 
J dbja(uj) — D t ot = Tre 2 n/m, where n is the density 
and m is the mass of the electrons, is exhausted entirely 
by the delta function Drude contribution alone. When 
the material becomes dirtier, the Drude peak of the opti- 
cal conductivity acquires the Lorentzian shape with the 
width of 1/r. The conductivity sum rule is still ex- 
hausted by the Drude part alone when only the impurity 
scatterings are present in the system. The total weight 
D to t, however, can change as the inpurity scatterings or 
other interactions are introduced when we consider a fi- 
nite bandwidth, because the projection to a restricted 
basis set disregards all excitations to higher energy than 
the bandwidth. When other interactions are present, the 
free-carrier Drude weight is reduced by the quasiparti- 



cle renormalization factor Z such that D = Z~ 1 D to t, 
and the missing spectral weight from the Drude part is 
transfered to a higher energy region of cr(u)) reflecting the 
excitation of incoherent scatterings. 

The experimentally measured <j(w) in the normal state 
A3C60 shows a remarkable reduction of Drude weight 
and, concomitantly, a pronounced MIR absorption below 
the inter-band absorption peak: DeGiorgi et al. found a 
pronounced MIR peak around 0.06 eV and analysized 
that the Drude weight is reduced to about 0.1 — 0.2 of 
the total intra-band spectral weight while Iwasa et al. 
observed the MIR absorption peak around 0.4 eV and de- 
termined that the Drude weight is reduced to about 0.6 of 
the total intra-band spectral weight Although their 
results show somewhat different Drude weight and MIR 
absorption energy each other, the pronounced suppres- 
sion of the Drude weight and the accompanying MIR ab- 
sorption imply strong electron-phonon and/or electron- 
electron interactions in this material. 

In order to understand this unusual feature in a(u)) 
of doped fullerenes, Gunnarsson et al. studied the ef- 
fects of the electron-phonon interaction on <j{ui) assum- 
ing that the Migdal theorem is valid ||. They showed 
that the electron-phonon interaction leads to a narrow- 
ing of the Drude peak by the factor Z = 1 + A, where A 
is the dimensionless electron-phonon coupling constant, 
and a transfer of the depleted Drude weight to a MIR re- 
gion at somewhat larger energies than the phonon energy. 
Their results, however, are far from sufficient to describe 
experimental observations. Therefore, they hinted that 
the Coulomb interaction between conduction electrons, 
which is neglected in their study, could lead to futher re- 
duction of the Drude weight and more pronounced MIR 
absorption. On the other hand, one of the present au- 
thors recently found, by studying the NMR coherence 
peak supression in the fullerene superconductors, that 
the Coulomb interaction between conduction electrons, 
characterized by UNp ~ 0.3 — 0.4, where U is the effec- 
tive Coulomb interaction and Np is the density of states 
(DOS) at the Fermi level, should be included in addition 
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to the electron-phonon interaction to understand the var- 
ious experimental observations in fullerenes in a coherent 
way We, therefore, included the electron-electron as 
well as electron-phonon interactions at the presence of 
the impurity scatterings in the present paper, to better 
understand the experimentally obsered unusual features 
in the optical spectra of the fullerene superconductors in 
the normal state. 

For fullerene superconductors, the Fermi energy ep = 
S/2 k 0.2 — 0.3 eV and the average phonon frequency 
uj pfl « 0.05 — 0.15 eV, where B is the bandwidth. There- 
fore, LUph/eF ~ 1 for fullerenes unlike conventional met- 
als, where ujph/tF ^ 1- When uj p h/sF ~ 1, the phonon 
vertex correction becomes important because the Migdal 
theorem does not hold Jl-Q, and the frequency de- 
pendence of the effective Coulomb interaction, V e ff(u>), 
should be considered because the frequency scale at 
which V e ff(w) varies is comparable with that of electron- 
phonon interaction In this present work, concerned 
with the effects of the Coulomb and electron-phonon 
interactions on the optical spectra in the narrow band 
fullerene superconductors, the vertex correction is incor- 
porated in calculating the electron self-energy P,p|,p^| . 
The Coulomb interaction, modelled in terms of the onsite 
Hubbard repulsion, is included on an equal footing with 
the electron-phonon interaction, and considered fully self- 
consistently in calculating the effective electron-electron 
interaction |§,|l0|. The effective electron-electron inter- 
action becomes frequency dependent through the screen- 
ing. The impurity effects are included with the t-matrix 
approximation. 

Through the relation 



H{ip) = G- 1 {ip)-G^\ip), 



(1) 



one obtain the electron self-energy T,(ip) in the Mas- 
tubara frequency, which gives in the real frequency 
after the analytic continuation. Go and G are, respec- 
tively, the bare and renormalized electron Green's func- 
tions. S(u;) or Z(lu), where the renormalization function 
Z(u>) is given by = u> — ujZ{uj), defines the single- 

particle Green's function of an interacting system as 



G 1 = uj - £, k - £(ci>) = ujZ(uj) - £ fc 



(2) 



where is the electron energy measured from the chemi- 
cal potential, £/. = Ek — Then, the optical conductivity 
can be obtained by calculating the current-current cor- 
relation function, TL(ioj), using the renormalized Green's 
function obtained from solving Eq. (Q) self-consistently. 
The calculated optical conductivity shows a strong re- 
duction of Drude weight and a broad MIR absorption, 
although the MIR feature around 0.06 eV is less pro- 
nounced and broader compared with experimental ob- 
servations. 

This paper is organized as follows: In the follow- 
ing section, we present the Eliashberg-type formalism 
in the Matsubara frequency to calculate the renormal- 
ized Green's function with the impurity, electron-phonon, 



and Coulomb interactions included self-consistently. We 
then describe the analytic continuation procedure to ob- 
tain the renormalization function Z(u>) in the real fre- 
quency. The optical conductivity calculated with the 
renormalized Green's function is presented in Sec. III. 
We will discuss how the Drude part and the MIR ab- 
sorption of <j(oj) are affected as the impurity scattering 
rate, the electron-phonon and electron-electron interac- 
tions are varied. These result will then be compare with 
the experimental observations. Finally, Sec. IV is for the 
summary and some concluding remarks. 



II. FORMALISM 

The optical conductivity is calculated from the 
current-current correlation function, H(u>), as a(uj) — 
~ lim 9 ^ H(q,uj) [jll). We use the approximation where 
the electron self energy is momentum independent. In 
this case, it can be shown that the vertex correction 
in the current-current correlation function vanishes for 
q -> Q. This leads to 

2 e 2 i 

n(iuj m ) = Q _ 2T/ G & iPn + iuJ m )G(j), ip n ), 



P W-n 



(3) 



where ip n = 7rT(2n + 1) and iuj m = 2itTm are, re- 
spectively, fermion and boson Mastubara frequencies, 
where T is the temperature, m and n are the integers. 
f3 = l/ksT, and V is the volume. The evaluation of Eq. 
(ft) using Eq. (0) produces 



2ne 2 n 1 v-^ 

nW = — 7. E 



8(ip + iuj) — 9{ip) 

m P (p + uj)Z{ip + iuj) — pZ{ip) 

ip 



(4) 



in the Mastubara frequency. After performing the ana- 
lytic continuation of itu — » uj + i5, to the real frequency, 
the optical conductivity is given by 



cr(w) 



1 e 2 n 



x Re 



de[f F (e)-f F (e + uj)} 

9{e+i5) ~0(e+uj+i5) 
eZ(e + iS) - (e+uj)Z(e+uj + iS) 
0(e-i5)-0(s+uj + i5) 



sZ(e-iS) - {e+uj)Z(e+uj+iS) 



where Jf{e) — l/(l + e' 3£ ) is the Fermi distribution func- 
tion, and 8(uj + iS) = tan -1 [ li)Z ^ f f ,^ ] ■ The finite con- 
duction bandwidth B with a constant DOS is explicitly 
considered through the factor of 6, which is ir/2 for the 
usual case of infinite bandwidth metal. In order to calcu- 
late the optical conductivity from Eq. (||) we need Z(uj) 
which defines single-particle interacting Green's function 
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G(uj). This can be obtained by solving Eq. ([j]) self- 
consistently. The electron self-energy is obtained by cal- 
culating the exchange diagram of the renormalized elec- 
tron Green's function and the effective electron-phonon 
and Coulomb interactions with the vertex correction in- 
cluded via the method of Nambu. The Coulomb interac- 
tion, modelled in terms of the onsite Hubbard repulsion 
for simplicity, is included on an equal footing with the 
electron-phonon interaction. The impurity effects are in- 
cluded with the i-matrix approximation. The Eliashberg- 
type equation can be written in the Mastubara frequency 
as 
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ZnPn = Pn + 7j / y [Vi.(" ~ m ) ~ Kh{n ~ m) 

m 

+ X sp (n - m)]29 m T + —9 n , 



(6) 



where 

kernel. X c h(n 



tan _1 ( 



and A 



ph{ 



m) 



a 2 F(n)2n 

[Q= + bn-Pm) 3 ] 



is the electron-phonon interaction 

and X sp (n — m) are, respectively, 
the interactions in the charge and spin channels due 
to the Hubbard repulsion. They are determined self- 
consistently as 

\ ch (k) = UN F i~-Xn+xlH± + VXn]\, (?) 



X sp (k) = UN F \-+ X n + xl Ml - 1/Xr. 



(8) 



where Xn(k) is the dimensionless susceptibility given by 



Xn(k) 



N F U 1 



e F (3 



'ivi+k- 



(9) 



The r on the right hand side of Eq. (||) represents the ver- 
tex correction satisfying the Ward- identity [jl3|. When 
we neglect the vertex correction, T = 1. If we assume 
a weak frequency dependence of T, the vertex function 
r reduces to Z(ip m ). In this work, we treat the vertex 
correction exactly, and T is given by 



ip n Z(ip n ) - ip m Z(ip n 



(10) 



Solving Eq. (g) self-consistently yields Z(iu>) in the Mas- 
tubara frequency. In order to calculate cr(cj), analytic 
continuation of iuj — > uj + i8 should be performed to get 
Z(uj) in real frequency. The numerically exact analytic 
continuation of standard Eliashberg equation is usually 
performed by the iterative method developed by Mar- 
siglio, Shossmann, and Carbotte (MSC) using a mixed- 
representation |]l4| . But when we include the vertex func- 
tion exactly, the MSC method can not be applied because 
it needs a specific form of equation. Here, in order to con- 
sider vertex correction exactly, we do the alanytic con- 
tinuation by employing the iterative method extended 



by Takada Jig] . In this case, Eq. (g) is transformed to a 
mixed representation as fallow: 

z(uj) = z{uj) + dnp(n)^[n B (n) + n F (uj + n)} 
-{oj+n)z(uj+n)-ujZ{ojh 



x G(uj + fl) 



+ [n B (ty+n F (tt-uj)]G(ijj-n) 
' (uj — Q)Z (lu — CI) — ujZ(uj) 



-Cl 



where 



1 f°° / 1 

z{u) = i+ ^2 dnp(n) (- 

\ip m Z{ip rn )-ujZ(uj) 
X G[ip m ) I- 



(11) 
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ip m - UJ 

P(CL) = —ImMCL) 

A(O) = x ch {ci) - x ph (ci) - x sp (ci) 

G{i Pm ) = 2e(ip m ), G r {uj) = 2i0[u). 



Cl ip m — uj + Cl / 
i 6{uj) 



(12) 



Z(uj) of Eq. ( |TT| ) represents the renormalization function 
obtained by substituting iuj to ui + id before the fre- 
quency summation. The second term is the correction to 
Z(uj) to yield the correct retarded renormalization func- 
tion Z(uj) one would have obtained if the analytic contin- 
uation were performed after the frequency summation. 
Putting the solutions of Eq. (^J), Z(iuj), into the Z(iuj) 
Eq. ( pT| ) yields a self-consistent Eliashberg-type equation 
in the real frequency. Then, Z(uj) can be obtained by 
computing iteratively Eq. ( |ll| ) . 

In order to model fullerene superconductors, three 
truncated-Lorentization functions were used to represent 
o?F{Cl) as follow |JlO): 



?F(n) = Y. a l F ^ n )> 

v=l 

f v {ci) = ( i[(n-^V 

I 0, other wi 



+r 2 rs+r 2 
se. 



for\Cl — uj v \ < r c 



where F v (Cl) is the truncated Lorentizian centered at 
to v with the width of T = T c is the cutoff fre- 

quency of r c = 3r, and R is normalization constant 
such that f Q dflF ll (Cl) — 1. Various theoretical and 
experimental estimates do not agree well each other in 
terms of distribution of coupling strength o? v among dif- 
ferent modes. These estimates show, however, that the 
phonon frequency derived from intramolecular A g and 
H g modes are distributed over 0.03 — 0.2 eV with the 
total A in the range of 0.5 — 1 eV. In view of this, we 
represent the phonon modes with three groups centered 
around uj v = 0.04, 0.09, 0.19 eV, and 2N F a 2 v /LO v = 
0.3A S , 0.2A S , 0.5A S , respectively, for v = 1,2,3. Note 
that J2v=i ^FO? v jijJ v = X s . The A s sets the strength 



(13) 
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of a 2 F(fl) and N F a 2 F(£l)/\ s is independent of A s . For 
infinite bandwidth superconductors, A is equal to X s in 
the limit T — > 0. For a finite bandwidth system, how- 
ever, A is reduced from X s because the available states 
to and from which quasiparticles can be sccattered are 
restricted as the bandwidth is reduced. 



III. RESULTS 

The self-consistent equation of Eq. (|l|) is solved nu- 
merically as described in the previous section to obtain 
Z(uj). Then, the optical conductivity is calculated from 
the Eq. (^). Fig. 1 shows the optical conductivity <j{uj) 
as A is varied when the Coulomb interaction U is set 
to for a reference. Here, the Fermi energy ef, tem- 
perature T and impurity scattering rate 1/r are set to 
0.25, 0.001, 0.01 eV, respectively. This result shows 
quite a similar behavior to van den Brink et al. calcula- 
tion. As A is increased, the width of Drude peak becomes 
narrower and it's weight is transferred to a mid-infrared 
spectrum. However, the reduction of Drude weight is less 
than the factor of (1 + A), because of the finite bandwidth. 
The inset shows a MIR absorption spectra obtained by 
extracting Drude part from the total optical conductivity. 
In determining the Drude weight, fitting procedure was 
carefully employed and confirmed by examining zero fre- 
quency extrapolation in the Mastubara frequency which 
is proposed by Scalapino et al fl6|| . The three Lorentizian 
peaks of a 2 F(uj) in the electron-phonon paring kernel are 
attributted to the development of these MIR peaks. But, 
the MIR peaks are broadened and move to slightly high 
frequencies. Fig. 2 shows the MIR absorption due to the 
Coulomb interaction. The MIR part is also extracted by 
fitting as shown in the inset. In order to focus on how U 
affects the total optical conductivity, A is set to 0. Ef, 
T and 1/r are same as in Fig. 1. The Coulomb interac- 
tion induces the strong to dependence of renormalization 
function Z(uj), and the low frequency strong to depen- 
dence of Z(ui) distorts the Drude part of optical con- 
ductivity and induces the MIR absorption in the fairly 
low frequency region. As the impurity effect is enhanced, 
the MIR absorption due to Coulomb interactions tends 
to shift to higher frequency and finally merge together 
with the MIR peaks developed by electron-phonon inter- 
action, as shown in Fig. 3 for UNp — 0.3 and A = 0.7. 
Note that the position of this merged MIR peak in Fig. 
3 is around and above 0.06 eV which is experimentally 
observed value of DeGiorgi. 

Fig. 4 is a(uj) of doped fullerenes with T = 0.005 eV, 
e f = 0.25 eV, 1/r = 0.1 eV, UN F = 0.3, and A = 0.7, 
which is to be compared with the experimental observa- 
tions. The Drude weight is reduced to 0.467 of the total 
intra-band optical weight. The reduction factor of the 
Drude weight by electron phonon interaction is 1 + A, 
and the finite bandwidth futhcr restricts the reduction 
factor. It therefore seems unlike that the Drude weight 



less than about 0.6 of the total intra-band spectral weight 
can be explained without the Coulomb interactions, when 
we take A ~ 0.7 — 0.8. The Coulomb interaction sup- 
presses the Drude part substantially by inducing u> de- 
pendence of the renormalization function Z(lu) in the low 
frequency region. We think that the large reduction of 
Drude weight like the DeGiorgi experiment is a result 
of the strong Coulomb interaction between conduction 
electrons in addtion to the electron-phonon interaction. 
However, our results are still not sufficient to explain ex- 
perimentally founded results: (a) The Drude weight is 
about 0.46 of the total intra-band optical weight with a 
resonable set of parameter values while DeGiorgi found 
0.1 — 0.2. (b) The MIR absorption is very broad which 
begins around 0.02 eV, has a peak around 0.07 eV and 
extends well over the Fermi energy. 



IV. SUMMARY AND CONCLUSION 

In this paper, we tried to give an explanation for the 
unusual behavior of optical conductivity in the normal 
state A^Ceo- It is generally accepted that the fullerene 
superconductor could be characterized by the phonon- 
mediated s-wave superconductor 0J^]. However, a few 
experiments like optical conductivity still remain not 
understood by the electron-phonon scattering together 
with the disorder effects. Our motivation lies in that 
the fullerene superconductors have such a narrow band- 
width that the phonon frequency, the Coulomb interac- 
tion, and the Fermi energy are all comparable, uj p h ~ 
V ~ E F - In order to consider properly the Coulomb 
interaction and the narrow bandwidth of fullerene super- 
conductors, the self-consistent Eliashberg-type coupled 
equations are solved to obtain the renormalized Green's 
function. The theory includes the frequency dependent 
screened Coulomb interaction together with the electron- 
phonon interaction and includes the vertex correction via 
Nambu's method. In order to treat the vertex function 
exactly, analytic continuation is performed via the it- 
erative method of mixed repersentation which is devel- 
oped by Takada. Once we get renormalization function 
Z(lo) in real frequency, we can calculate optical conduc- 
tivity in normal states. As we expected, the electron- 
phonon interaction is not suffficient to resolve the sub- 
stantial reduction of Drude weight and pronounced MIR 
peak. The strong Coulomb interaction induces to depen- 
dence in renormalization function Z(lu). As a result, the 
Drude form in optical spectra is distorted accompanying 
the reduction of Drude weight. When the impurity ef- 
fect is enhanced, the MIR absopption induced by strong 
Coulomb interaction merge together with the MIR peaks 
due to electron-phonon scattering showing large reduc- 
tion of Drude weight and MIR peak around 0.06 eV. 
Although it is not sufficent to explain experimentally 
founded results, our result is close to DeGiorgi's experi- 
ment. We improve Gunnarsonn's calculation by consid- 
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ering the electron-electron interaction and finite band- 
width effects explicitly. In conclusion, the unusual be- 
havior of optical conductivity of the normal state A3C60 
reveals the fact that both the Coulomb interaction and 
electron-phonon interaction are important in examining 
dynamical properties of fullcrcne superconductors. 
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Figure Captions 



Figure 1. The optical conductivity as a function of w 
for various electron-phonon coupling constants A when 
T = 0.001 eV, e F = 0.25 eV and 1/r = 0.01 eV. U is set 
to for a referrence. As A is increased, the width and the 
weight of the Drude peak are reduced. The inset shows 
a decomposition of the total conductivity into the Drude 
and MIR parts for A = 0.7. 



Figure 2. The MIR spectra induced by the Coulomb 
interaction when T = 0.001 eV, e F = 0.25 eV, 1/r = 0.01 
eV and A = 0. The inset shows a decomposition of the 
total conductivity, as in Fig. 1, into the Drude and MIR 
parts for UN F = 0.5. 

Figure 3. The MIR spectra as the impurity scattering 
rates 1/r are varied when T = 0.005 eV, e F — 0.25 eV, 
A = 0.7 and UN F = 0.3. When 1/r = 0.01 eV the lower 
peak is mainly from the Coulomb interaction while the 
other peaks are from the electron-phonon interaction. As 
1/r is increased, these peaks are merged altogether and 
finally evolve into a single broad peak around 0.06 ~ 0.1 
eV. 

Figure 4. The total optical conductivity with it's 
Drude and MIR parts for T = 0.005 eV, e F = 0.25 eV, 
1/r = 0.1 cV, A = 0.7 and UN F = 0.3. The Drude part 
in the low frequency region is substantially suppressed 
due to the strong Coulomb interaction. Consequently, 
the missing spectral weight is transfcrcd to the broad 
MIR peak, which peaks around 0.07 eV and extends well 
into the higher energy region. The ratio of the MIR 
spectral weight to the total intra-band spectral weight 
is 0.533. 
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